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G\ ' Abstract 



. We reconsider the time-dependent Born-Oppenheimer theory with the goal to carefully 

separate between the adiabatic decoupling of a given group of energy bands from their 
orthogonal subspace and the semiclassics within the energy bands. Band crossings are allowed 
and our results are local in the sense that they hold up to the first time when a band crossing 
is encountered. The adiabatic decoupling leads to an effective Schrodinger equation for the 
f^S nuclei, including contributions from the Berry connection. 



"?h 1 1 Introduction 
pk 

Molecules consist of light electrons, mass m e , and heavy nuclei, mass M which depends on 
the type of nucleus. Born and Oppenheimer || wanted to explain some general features of 
molecular spectra and realized that, since the ratio m e /M is small, it could be used as an 
expansion parameter for the energy levels of the molecular Hamiltonian. The time- independent 
Born-Oppenheimer theory has been put on firm mathematical grounds by Combes, Duclos, and 
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^ Seiler || , Hagedorn ||| , and more recently in |l6| 



With the development of tailored state preparation and ultra precise time resolution there is 
a growing interest in understanding and controlling the dynamics of molecules, which requires 
an analysis of the solutions to the time- dependent Schrodinger equation, again exploiting that 
m e /M is small. The molecular Hamiltonian is of the form 

H = 2^" ( " iVx ~ Aext(x) ) 2 + m ( " iVx + Acxt 2 + Vc ^ + v ™( x > x ) + K w • w 

For notational simplicity we ignore spin degrees of freedom and assume that all nuclei have 
the same mass. We have k electrons with positions {x\, . . . , Xk} = x and I nuclei with positions 
{Xi, . . . , Xi} = X . The first and second term of H are the kinetic energies of the electrons and of 
the nuclei, respectively. An external magnetic field is included through the vector potential A ex i. 
Electrons and nuclei interact via the static Coulomb potential. Therefore V e is the electronic, 
V n the nucleonic repulsion, and V en the attraction between electrons and nuclei. V e and V n may 
also contain an external electrostatic potential. 
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Figure 1: The schematic spectrum of H e (R) for a diatomic molecule as a function of the sepa- 
ration R of the two nuclei. 

In atomic units (m e = h = 1) the Hamiltonian (|l|) can be written more concisely as 



emphasizing that the nuclear kinetic energy will be treated as a "small perturbation". H C (X) 
is the electronic Hamiltonian for given position X of the nuclei, 



H C (X) is a self-adjoint operator on the electronic Hilbert space L 2 (M 3fc ) restricted to its antisym- 
metric subspace. Later on we will need some smoothness of H e (X), which can be established 
easily if the electrons are treated as point-like and the nuclei have an extended, rigid charge 
distribution. 

Generically H e (X) has, possibly degenerate, eigenvalues E\(X) < E2(X) < . . . which termi- 
nate at the continuum edge Thereby one obtains the band structure as plotted schemat- 
ically in Figure 1. The discrete bands Ej(X) may cross and possibly merge into the continuous 
spectrum as indicated in Figure 2. 

Comparing kinetic energies, we find for the speeds \v n \ ~ (m c /Af) 1//2 |u c |, which means that 
on the atomic scale the nuclei move very slowly. If we regard X(t) as a given nucleonic trajectory, 
then H e (X(t)) is a Hamiltonian with slow time variation and the time-adiabatic theorem [fi5| , |l4| , 
|J can be applied Q. For us X are quantum mechanical degrees of freedom. The Hamiltonian H 
of (||) is time-independent and we can only exploit that the nucleonic Laplacian carries a small 
prefactor. To distinguish, we refer to our situation as space-adiabatic. Since the nuclei move 
very slowly, their dynamics must be followed over sufficiently long times. From the speed ratio 
we conclude that these times are of order (m e /M) l l 2 in atomic units. To simplify notation we 




iV x - ^ext(x)) + V e (x) + V en (X, x) + V n (X) . 



(3) 



define 




(4) 



as the small dimensionless parameter. Then 




(5) 



2 



and we want to study the solutions of the time-dependent Schrodinger equation 



ie^ = (6) 

in the limit of small e. 

The crude physical picture underlying the analysis of (||) is that the nuclei behave semiclas- 
sically because of their large mass and that the electrons rapidly adjust to the slow nucleonic 
motion. Thus, in fact, the time-dependent Born-Oppenheimer approximation involves two lim- 
its. If the electrons are initially in the eigenstate Xj(Xo) of the j-th band with energy Ej(Xo), 
where Xq is the approximate initial configuration of the nuclei, then the j-th band is adiabati- 
cally protected provided there is an energy gap separating it from the rest of the spectrum. Thus 
at later times, up to small error, the electronic wave function is still in the subspace correspond- 
ing to the j-th band. But this implies that the nuclei are governed by the Born-Oppenheimer 
Hamiltonian 

H £ bo = <?\ ( - + A cxt (X)) 2 + Ej{X) . (7) 

Since e< 1, ^bo can ^ e anar yzed through semiclassical methods where to leading order the 
contributions come from the classical flow $* corresponding to the classical Hamiltonian H^ = 
t;P 2 + Ej(q) on nucleonic phase space. 

In general, Ej(X) may touch another band as X varies. To allow for such band crossings we 
introduce the region A C M n , n = 31, in nucleonic configuration space, such that Ej restricted 
to A does not cross or touch any other energy band. The classical flow $* then has A x R n as 
phase space and is defined only up to the time when it first hits the boundary dA x W 1 . Up to 
that time (0) still correctly describes the quantum evolution. To follow the tunneling through 
a band crossing other methods have to be used [[ll], M, in particular, the codimension of the 
crossing is of relevance. 

The mathematical investigation of the time-dependent Born-Oppenheimer theory was ini- 
tiated and carried out in great detail by Hagedorn. In his pioneering work || he constructs 
approximate solutions to (g) of the form 4> q (t), P (t) ® Xj(,<l(t))i where 4> q (t), P (t) ls a coherent state 
carried along the classical flow, (q(t),p(t)) = <3?*(go,R>)- The difference to the true solution 
with the same initial condition is of order ^fe in the L 2 -norm over times of order e _1 in atomic 
units and the approximation holds until the first hitting time of dA x R n . In a recent work 
Hagedorn and Joye [10] construct solutions to (@) satisfying exponentially small error estimates. 



In Hagedorn 's approach the "adiabatic and semiclassical limits are being taken simultaneously, 



and they are coupled [10]". 

In our paper we carefully separate the space- adiabatic and the semiclassical limit. One 
immediate benefit is the generalization of the first order analysis of Hagedorn from coherent 
states to arbitrary wave functions. 

Let us explain our result for the space-adiabatic part in more detail. We assume that there 
is some region A C W 1 in the nucleonic configuration space, such that some subset cr*(X) of 
a(H c (X)) is separated from the remainder of the spectrum by a gap for all X £ A, i.e. 

dist(cr*(X), a(H c (X)) \ a*(X)) > d > for all X G A . 

A could be punctured by small balls (for n = 2) because of band crossings. A could also terminate 
because the point spectrum merges in the continuum, which physically means that the molecule 
loses an electron through ionization. Let P*(X) be the spectral projection of H C (X) associated 
with <t*(X) and = dX P*{X). We will establish that the unitary time evolution e ~ lHet / s 
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Figure 2: The wave function can leave RanP* in two different ways. Either by transitions to 
other bands (a) or through the boundary of A (b). 



agrees on RanF* with the diagonal evolution e _J ^ ia s*^ e generated by #J iag := P*H e P* up to 
errors of order e as long as the leaking through the boundary of A is sufficiently small. 

To complete the analysis one has to control the flow of the wave function through <9A. One 
possibility is to simply avoid the problem by assuming that A = R n , hence <9A = 0. We will 
refer to this case as a globally isolated band. Of course, the set {(X,y) £ R n x R : y £ a*(X)} 
may contain arbitrary band crossings. As one of our main results, we prove that the subspace 
Rani 3 * is adiabatically protected. In particular for the purpose of studying band crossings the 
full molecular Hamiltonian may be replaced by a simplified model with two bands only. 

In general one has dA ^ 0, to which we refer as a locally isolated band. To estimate the 
flow out of A the only technique available seems to be semiclassical analysis. But this requires a 
control over the semiclassical evolution, for which one needs, at present, that {(X,y) 6 AxM: 
y £ a*(X)} contains no band crossings. Then {(X, y) £ A x R : y £ a^(X)} = Uj{(X,y) £ 
A x R : y = Ej(X)} is the disjoint union of possibly degenerate energy bands Ej(X). We will 
prove that each band separately is adiabatically protected. 

In the special case where a*(X) = Ej(X) is a nondegenerate eigenvalue for X £ A, e~ tHdia ^ t ^ e 
is well approximated through e _J ^So*/ £ Q n L 2 (R n ). Since H^ is a standard semiclassical 
operator, one can easily control the X-support of the wave function and therefore prove a result 
for rather general A C R n , for the details see Theorem Roughly speaking, it says that if (j)t is 
a solution of the effective Schrodinger equation for the nuclei 

is^ = H^ t , (8) 
with supp^o C A, then, modulo an error of order s, 

is a solution of the full Schrodinger equation (||) with initial condition ipo(X, x) = <po(X)xj{X, x) 
as long as 4>t is supported in A up to L 2 -mass of order e. This maximal time span can be 
computed using the classical flow <!>*. 

As first observed by Mead and Truhlar , in general H^ Q acquires as a first order correction 
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an additional vector potential A geo (X) = — i(xj(X), VxXj(X)) and (|7|) has to be replaced by 

1 



Hl i0 = e z -i K -iVx + A cxt (X)+A gC0 (X)) + Ej(X) . (9) 

Multiplying Xj(X) with a smooth X-dependent phase factor induces a gauge transformation 
for ^geo, which implies that the physical predictions based on (||) do not change, as it should 
be. As noticed in JO]], if A is not contractible, then ^4g eo cannot be removed through a gauge 
transformation and (0) and (0) describe different physics. Berry realized that geometric phases 
appear whenever the Hamiltonian has slowly changing parameters. Therefore A geo (X) is referred 



to as Berry connection, cf. [22] for an instructive collection of reprints. In fact, the motion 
of nuclei as governed by the Born-Oppenheimer Hamiltonian (g) is one of the paradigmatic 
examples for geometric phases. 

If <7*(A) = E(X) is A;- fold degenerate, not much of the above analysis changes. H^ becomes 
matrix-valued and acts on L 2 (M. n )® k , i.e. 

H £ bo= (j { - iVx + A cxt (X)Y + EjiX)^ i kxk+ £ -^-ieV x yA geo (X)+A gco (Xy(-ieVx) 

The connection A geo (X) contains in general also off-diagonal terms and matrix- valued semiclas- 
sics must be applied. However, since the only nondiagonal term is in the subprincipal symbol, 
the leading order semiclassical analysis reduces to the scalar case and, in particular, agrees with 
the nondegenerate band case. We do not carry out the straightforward extension of Theorem 
H below to the degenerate band case, because the technicalities of matrix-valued semiclassics 
would obscure the simple ideas behind our analysis. 

In their recent work Martinez and Sordoni independently study the time-dependent 
Born-Oppenheimer approximation as based on techniques developed by Nenciu and Sordoni 
pel . They consider the case of a globally isolated band for a Hamiltonian of the form (|l|) with 
smooth V and A ext = 0. They succeed in proving the adiabatic decoupling to any order in e for 
subspaces which are e-close to the unperturbed subspaces P* considered by us. With this 
result, in principle, higher order corrections to the effective Hamiltonian ([?]) could be computed. 

The paper is organized as follows. Section 2 contains the precise formulation of the results. 
Section 3 gives a short discussion of the semiclassical limit of H^ and on how such results 
extend to the full molecular system. Proofs are provided in Section 4. In spirit they rely on 
techniques developed in |23|| in the context of the semiclassical limit for dressed electron states. 
In practice the Born-Oppenheimer approximation requires several novel constructions, since the 

2 

"perturbation" — ^-A increases quadratically. 

Our results can be formulated and proved in a more general framework dealing with, possibly 
time-dependent, perturbations of fibered operators. Also the gap condition can be removed by 
using arguments similar to those developed by Avron and Elgart in |TJ. The general operator 



theoretical results will appear elsewhere [24| 



2 Main results 

The specific form (||) of the electronic part of the Hamiltonian will be of no importance in the 
following. Thus we only assume that 

H e = 1^ dX H C (X) , H C (X) = H c0 + H cl (X) , 



5 



where H e Q is self-adjoint on some dense domain D C TC e and bounded from below and H C \(X) G 
C(TC e ) is a continuous family of self-adjoint operators, bounded uniformly for X G W 1 . Thus 
F e is self-adjoint on D(H e ) = L 2 (R n ) ® V C H := L 2 (R n ) ® W e and bounded from below. 
For the definition of L 2 (R n ) we equip T> with the graph-norm || • ||# c0 , i.e., for ip £ T>, 
U\\ Hc0 = \\H e0 i;\\ + U\\. 

Let A cxt G C£(R n ,K n ), where for any open set O C R m , m G N, C§(fi) denotes the set of 
functions / G C fc (0) such that for each multi-index a with \a\ < there exists a C a < oo with 

sup\d a f(x)\ <C a . 

Then ^ ( - iV x + AsctPO) is self-adjoint on M^ 2 (lR n ), the second Sobolev space, since — iV x 
is infinitesimally operator bounded with respect to —Ax- It follows that 

# £ = y(-iVx + A cxt {X))\l + H c (10) 

self-adjoint on D{H e ) = W 2 {R n ) ®H c n D(H C ). 

For X G A, A C M n open, we require in addition some regularity for H C (X) as a function of 

X: 

H k flel(-)eCb(A,^(We)). 

The exact value of will depend on whether A = R n or A C M n . For the type of Hamiltonian 
considered in the introduction, cf. (|l]), all the above conditions including Condition are 
easily checked and put constraints only on the smoothness of the external potentials and on the 
smoothness and the decay of the charge distribution of the nuclei. For point nuclei fails and 
a suitable substitute would require a generalization of the Hunziker distortion method of jjHJ . 

We will be interested in subsets of {(X, s) G A x R : s G cr(H c (X)} which are isolated from 
the rest of the spectrum in the following sense. 

S For X G A, let a*(X) C a{H c {X)) be such that there are functions f± G Cb(A, M.) and a 
constant d > with 

[/_(*) + d, f + (X) -d]H a*(X) = a*(X) 

and 

[/_(*), f+(X)] n (a(H c (X) \ a m (X)) = . 

We set P* = dX P*(X), where P*(X) = t af rx)(H e (X)) is the spectral projection of 
H C (X) with respect to cr*(X). As explained in the introduction we have to distinguish two 
cases. 

(i) Globally isolated bands 
We assume A = R n and let 

tf| iag := P* H £ P, + H £ . (11) 

Since we aim at a uniform result for the adiabatic theorem, we introduce the Sobolev spaces 
\y 1 ' s (W l ) and W 2 ' e (R n ) with respect to the e-scaled gradient, i.e. 

W 1,e (R n ) := {<P G L 2 (R n ) : U\\ m „ := ||e | || + \\cj>\\ < oo} 
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and 

W 2 ^(R n ) := {</> G L 2 (M") : ||0|| wa ,. := ||e 2 A0|| + \\4>\\ < 00} . 

Alternatively we will project on finite total energies and define £(H £ ) := l(_ 00 ,f](-^' e ) as the 
projection on total energies smaller than £. 

Theorem 1. Assume H3 and S /or A = W 1 . Then H^ iag is self-adjoint on the domain of H e . 
There are constants C,C < 00 such that for allt&M. 



e -iHH/e _ e -iH- diag t/s 



<eC(l + \t\) 3 (12) 

C(W 2 *®Hc,H) 



and for all £ EM. 

J (e-wn/e _ g-ifl^t/^ £ (H £ )j^ H < eC(l + \£\)(l + \t\). (13) 

C(W 2 ' 6 <8> TL e ,Tl) denotes the space of bounded linear operators from W 2,e (8> 7i e H equipped 
with the operator norm. 

This result should be understood as an adiabatic theorem for the subspaces Rani 3 * and 
RanP^, which are not spectral subspaces. Let us point out one immediate application of Theo- 
rem [|. The behavior near band crossings is usually investigated using simplified models involving 
only two energy bands and ignoring the rest of the spectrum, cf. 11, [7]]. Theorem [l] shows that 
this strategy is indeed justified modulo errors of order e. 

(ii) Locally isolated bands 

<J*{X) = E(X) is a nondegenerate eigenvalue for all X £ A. A may now be any open subset 
of W 1 and for such a A we assume and S. We also assume that A is connected. Otherwise 
one could treat each connected component separately. 

It is easy to see that, given Hqq and S, the family of projections P*(-) € Cg°(A, C(Ti e )). 
However, in order to "map" the dynamics from Rani 3 * to L 2 (A) we need in addition a smooth 
version x(') £ Cb°(^-)^e) of the normalized eigenvector of H e {X) with eigenvalue E(X). In 
other words we require the complex line bundle over A defined by P* to be trivial. This always 
holds for contractible A, but, as discussed below, also for some relevant examples where A is not 
contractible. 

Given a smooth version of x(X) with ||x(A)|| = 1, one has Re(x(X) ,V xx(X)) = 0, but, 
in general, lm{x(X), V xx{X)) 7^ 0. In the following we distinguish two cases: Either it is 
possible to achieve Im{x(X) ,V xx(X)) = by a smooth gauge transformation x(X) — ► x(X) = 
e^^xiX) or not. In the latter case 

A geo (X) :=-i( X (X),Vxx(X)) 

is the gauge potential of a connection on the trivial complex line bundle over A, the Berry 
connection, and has to be taken into account in the definition of the effective operator 

#bo == j(-iVx + A ext {X) + A geo {X)Y + E(X) (14) 

with domain W 2 (R n ). Thus A geo acts as an additional external magnetic vector potential. 
Although ^4 ext and A geo appear in H^ with an e in front only, and therefore are not retained 
in the semiclassical limit to leading order, they do contribute to the solution of the Schrodinger 
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equation for times of order e _1 . If the full Hamiltonian is real in position representation, as it 
is the case for the Hamiltonians considered in the introduction whenever A ext = 0, then x(X) 
can be chosen real-valued. If, in addition, A is contractible, the existence of a smooth version 
of x(X) with lm{x(X), Vxx(X)) = follows. 

To define #g on L 2 (R n ) through ([lj), the functions E(X) and A geo (X), which are a priori 
defined on A only, must be continued to functions on R n . Hence we arbitrarily extend E{X) 
and A geo (X) to functions in C£°(IR n ) by modifying them, if necessary, on A \ (A — 5/5) (cf. (|l7l) ) 
for some 5 > 0. The parameter 5 will be fixed in the formulation of Theorem ||| and will appear 
in several places. It controls how close the states are allowed to come to dA. 

The generic example for the Berry phase is a band crossings of codimension 2 (cf. [p^ , 
[TT| , [7]]). If E(X) is an isolated energy band except for a codimension 2 crossing, then A = 
IR n \ {closed neighborhood of the crossing} is no longer contractible, but the line bundle is still 
trivial. Although the underlying Hamiltonian is real, the Berry connection cannot be gauged 
away. Within the time-independent Born-Oppenheimer approximation Herrin and Howland [|l^] 
study a model with a nontrivial eigenvector bundle. 

With the fixed choice for x(X) we have 

RanP* = Ij® dX (f>(X)x(X); <j) £ L 2 (A) j C Tt . (15) 

Thus there is a natural identification IA : RanP* — > L 2 (R n ) connecting the relevant subspace 
on which the full quantum evolution takes place and the Hilbert space L 2 (R") on which the 
effective Born-Oppenheimer evolution is defined. According to (jig), we set 



U{cj)x) = ct>, i.e. (UP^)(X) = (x(X), (P^)(X)) Hc . 
Its adjoint U* : L 2 (R n ) — > RanP* is given by 



U*cf>= / dX HX)x(X) . 
Ja 

Clearly U is an isometry and U*U = 1 on RanP*. But U is not surjective and thus not unitary. 

By construction, e~ tH so t / e is a good approximation to the true dynamics only as long as the 
wave function of the nuclei is supported in A modulo errors of order e. Since is a standard 
semiclassical operator, the X-support of solutions of (||) can be calculated approximately from 
the classical dynamics generated by its principal symbol H c \(q,p) = ^p 2 + E(q) on phase space 
Z := R n x R n , 

Jt q=P ' Jt P= ~ VE ^- (16) 
The solution flow to (16) exists for all times and will be denoted by 



In order to make these notions more precise, we need to introduce some notation. The Weyl 
quantization of a £ C£°(Z) is the linear operator 

(a W ' £ </») (X) = (2ir)- n [ dYdka (?^-,e k] e -^ x - Y > k <P(Y) , 

as acting on Schwartz functions. a w,e extends to C(L 2 (W 1 )) with operator norm bounded 
uniformly in e (cf., e.g., Theorem 7.11 in ||). The wave functions with phase space support in 
a compact set r C Z do not form a closed subspace of L 2 (M n ). Hence we cannot project on this 
set. In order to define approximate projections, let for T C M. m , m S N, and for a > 



a := < z S T : inf \w — z\ > a > . (17) 

1 «>eiR m \r ' 
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Definition 2. An approximate characteristic function l(r» £ C£°(]R m ) of a set V C W 71 with 
margin a is defined by the requirement that l(r,o:) |r— a = 1 and U(r,«) |]R m \r = 0. 

If l(r,a) *s an approximate characteristic function on phase space Z , then the corresponding 
approximate projection is defined as its Weyl quantization l, r '^ . We will say that functions in 
RanlfJ 6 , have phase space support in T . 

For r C Z we will use the abbreviations 

T q := {q G 1" : (q,p) G T for some p G M™} , 
T p := {p G M n : (q,p) G T for some q G R"} . 

Let the phase space support T of the initial wave function be such that T q C A — 5. Then 
the maximal time interval for which the X-support of the wave function of the nuclei stays in 
A up to errors of order e can be written as 

/Lx(I\A) := [T 5 _(T,A),Tl(T,A)}, 

where the "first hitting times" T± are defined by the classical dynamics through 

I+(T,A) :=sup{i>0: ($ s (r)) g CA-i Vs G [0,t}\ 

and T_(r, A) analogously for negative times. This are just the first times for a particle starting 
in r to hit the boundary of A — 5 when dragged along the classical flow . 

The following proposition, which is an immediate consequence of Egorov's Theorem |4|, 
shows that for times in /^(r, A) the support of the wave function of the nuclei stays indeed 
in A — 5, up to errors of order e uniformly on RanlL r '^\ for any approximate projection l(p'^v 

Proposition 3. Let T C Z be such that T q C A — 5 and let ]L\_5 denote multiplication with 
the characteristic function of A — 5 on L 2 (R n ). For any approximate projection ^rp'a) an d an y 
bounded interval I C /^ ax (r, A) there is a constant C < oo such that for all t G I 



(1 - l A _i) e~ iH Bo^ l w > £ 



(r,«) 



< Ce. 

£(L 2 (R™)) 



W l {AiS) t^ £ a) UP*, where 1^ 



An approximate projection on T in Ji is defined as Pp 
is an approximate projection on T according to Definition [2] and 1(a,<5) is an approximate char- 
acteristic function for A. Using the latter instead of the sharp cutoff from IA* 
bounded set in W 2 ' £ (g> H e 



makes RanPp a 



whenever T p is a bounded set. 



Theorem 4. Assume Hqq and S with dim(RanP*(X)) = 1 for some open A C W 1 . Let T C Z 
be such that T q C A — 5 for some 5 > and T p bounded. For any approximate projection Pp 
and any bounded interval I C /^^(r, A) there is a constant C < oo such that for all t G I 



-iHH/e 



iH BO t / £ U 



< Ce. 



(18) 



Theorem ^ establishes that the electrons adiabatically follow the motion of the nuclei up to 
errors of order e as long as the leaking through the boundary of A is small. The semiclassics 
was used only to control such a leaking uniformly. However, for -fffj the limit e — > is a 
semiclassical limit and, as discussed in the following section, beyond the mere support of the 
wave function more detailed information is available. 
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3 Semiclassics for a single band 



The semiclassical limit of Equation (||) with a Hamiltonian of the form (14) is well understood 
and there is a variety of different approaches. For example one can construct approximate 
solutions ^> q {t) of (|8|) which are localized along a classical trajectory q(t), i.e. along a solution of 
(|l~6|) . Then it follows from Theorem |I| that 4> q (t)X is a solution of the full Schrodinger equation, 
@, up to an error of order e as long as q(t) £ A — 5. Roughly speaking, this coincides with the 
result of Hagedorn |J . In applications the assumption that the wave function of the nuclei is well 
described by a coherent state seems to be rather restrictive and a more general approach to the 
semiclassical analysis of a Schrodinger equation of the form @ is to consider the distributions of 
semiclassical observables, i.e. of operators obtained as Weyl quantization a W ' e of classical phase 
space functions a : Z —* R. 

Consider a general initial wave function (ft 6 £ L 2 (M n ), such that 4> e corresponds to a proba- 
bility measure p c \{dqdp) on phase space in the sense that for all semiclassical observables with 
symbols a e C£°(Z) 



lim 



,W,£ 



/ a(q,p) p c \(dqdp) 
Jz 



0. 



(19) 



The definition is equivalent to saying that the Wigner transform of (fi £ converges to p c \ weakly 



on test functions in C^(Z) [17]. An immediate application of Egorov's theorem yields 



lim 

£^0 



ls JH*t/e n W,e -iH^t/e ie\ 



(ao$*)(g,p) Pcl (dqdp) 







(20) 



uniformly on bounded intervals in time, where we recall that $* is the flow generated by (|l6|). 
In (p0|) one can of course shift the time evolution from the observables to the states on both 
sides and write instead 



lim 



,W,6 



a(q,p) Pd(dqdp,t) 



0. 



(21) 



Here (f>f = e~' lH ^ Bot ' £ (ff and p c \(dqdp,t) = (p c \ o ^ t )(dq dp) is the initial distribution p c \(dqdp) 
transported along the classical flow. Thus with respect to certain type of experiments the system 
described by the wave function <p £ behaves like a classical system. 

For a molecular system the object of real interest is the left hand side of pi] ) with <\>\ replaced 
by the solution ipf of the full Schrodinger equation and a w,e =: ag as acting on L 2 (M n ) replaced 
by a W ' e ®1 as acting on 7i. In order to compare the expectations of ag G with the expectations 
of a w,£ (g) 1, we need the following proposition. 



Proposition 5. In addition to the assumptions of Theorem let a S C^(Z) with 

d£ sup |£| |a (2) (x,£)| < do, 



/ 



(22) 



where denotes Fourier transformation in the second argument. Then there is a constant 

C < oo such that 



a 



W,e 



1 - U*a w ' £ U) t A - S P* || < Ce. 



For the proof of Proposition |5| see the end of Section 4.2. With its help we obtain the 
semiclassical limit for the nuclei as governed by the full Hamiltonian. 
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Corollary 6. Let F and I be as in Theorem 0. Let tp £ £ TL such that is satisfied for 
cf) e := UP*ip £ for some p c \ with suppp c \ C T — a. Let xfjf = e -iH e t/e^pE then for all a G C£°(Z) 
which satisfy Mi 



lim 



(f/4, (a W ' £ ®l)^ £ )- / a(q,p)p cl (dqdp,t) 



z 



(23) 



uniformly for t £ I. 



Translated to the language of Wigner measures Corollary |6| states the following. Let us 
define the marginal Wigner transform for the nuclei as 



W^MMp) ■■= (2^ n / dXe iX -v(i>?(q + eX/2),rt(q-£X/2)) 



We 



Then, whenever W £ nc (P*ipo)(q,p) dq dp converges weakly to some probability measure p c \(dq dp), 
^nuc(-f*V 't){l->P) dqdp converges weakly to (p c i o ^~ t )(dq dp). 

Corollary ^ follows by applying first Proposition || and then Theorem ^| to the left hand side 
in the difference (|23"|), where we note that lim e ^o ||(1 — P^ )ip £ \\ = an d thus also lim £ ^o ||(1 — 



Ph-S'^tW = f° r an y < This yields the left hand side of (pOj) and thus (23). 

We mention some standard examples of initial wave functions (jf of the nuclei which approx- 
imate certain classical distributions. The initial wave function for the full system is, as before, 
recovered as ip e = U*4> £ = <j) £ (X)x(X). In these examples one regains some control on the rate 



of convergence with respect to e which was lost in (19). 

(i) Wave packets tracking a classical trajectory. 
For 6 G L 2 (R n ) let 



Then \(j) £ ^ po (X)\ 2 is sharply peaked at go f° r £ small and its e-scaled Fourier transform is 
sharply peaked at pq. Thus one expects that the corresponding classical distribution is given by 
5(q — qo)5(p — po) dqdp. As was shown, e.g. in 23], this is indeed true for G L 2 (W 1 ) such that 



\p\4> G L 1 (R n ). Then Corollary | holds with (|D replaced by 

|(#, (a w > £ ®l)^ £ t )-a(q(t),p(t))\ 
= 0(y/e) (U\\ 2 L2 + U\\ L i \\\p\4\\ L i + \\\x\<t>\\ Ll U\\ L1 ) , (24) 



where (q(t),p(t)) is the solution of the classical dynamics with initial condition (qo,Po)- ( |24l ) 
generalizes Hagedorn's first order result in || to a larger class of localized wave functions. 

(ii) Either sharp momentum or sharp position. 
For (j) G L 2 (R n ) let 

<t>p (p) = 4>(p-Po/£) , 

where ^denotes the e-scaled Fourier transformation, then the corresponding classical distribution 
is p c i(dqdp) = 5(p — Po)\<fi(q)\ 2 dqdp. Note that the absolute value of <f> does not depend on e in 
that case. Equivalently one defines 
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and obtains p c \(dq dp) = 5(q — %)\<j){p)\ dq dp. In both cases one finds that the difference in (|2 
is bounded a constant times either e(||0[|^ 2 + \\4>\\i,i \\ IpI^Wl 1 ) f° r 0p o or e (ll ( / , ll|2 + II M<^||i 



for cf> £ qo . 

(Hi) WKB wave functions. 

For / G L 2 (R n ) and S G C l {R n ) both real valued let 

c/>%X) = f(X)e i ^ 1 , 

then p c \(dqdp) = f 2 (q) 5(p — VS(q)) dqdp. In this case one expects that (|2^ ) is bounded as yfe, 
which has been shown in |23(1 for a smaller set of test functions. 



4 Proofs 

4.1 Globally isolated bands 

We collect some immediate consequences of H3 and S. Using the Riesz formula 

P*(X) = ~ I dXR x (H e (X)) , (25) 

with 'y(X) a smooth curve in the complex plain circling <7*(X) only and R\(H C (X)) = (H C (X) — 
A) , one easily shows that -P*(-) G C^(M n , C{TL C )). Assumption S enters at this point, since it 
allows to chose j(X) locally independent of X. Hence, when taking derivatives with respect to 
X in (|25|), one only needs to differentiate the integrand. In particular one finds that 

P^{X){V X P*){X)P,{X) = 

/ d\R x (H c (X))P, ± (X)(V x H c )(X)R x (H c (X))P,(X). (26) 
2?r« J^(x) 



Since P*(X)(V X P*)(X)P*(X) = P^(X)(V X P*){X)P^(X) = 0, which follows from (V X P*)(X) 
(V X P?)(X) = (VxP*)(X)P*(X) + P.(X)(y x P*)(X), we have that 

(V X P*)(X) = PHX)(VxP*)(X)P,(X) + adjoint . (27) 



In (27) and in the following "+ adjoint" means that the adjoint operator of the first term in a 
sum is added. 

Starting with (|l2]), we find, at the moment formally, that 



,t/ £ 



i e 



f £ dse lH ^ s (H £ -H E diag )e- lH£s , (28) 
Jo 



where 

LIS L 

diag 



H £ ~ #dia e = P* H £ P* + adjoint 



P A 



y(-iVx+^extPO) 2 ,P* 



P* + adjoint . (29) 
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Let Da '■= -jV^ + A ext (X). Then the commutator is easily calculated as 
-2 i 

-ie(V x P*) • (eD A tg> 1) + 0{e 2 } 



(30) 

= -e(V x P*) ■ (eV x ®l) + 0(e 2 ), (31) 

where 0(e 2 ) holds in the norm of £(H, 7i) as e — > 0. For (|30| ) and ( |3l|) it was used that A ex t(X) 
and P*(X) are both differentiable with bounded derivatives and that A^^X) commutes with 
P*- 

Before we can continue, we need to justify (^|) by showing that is self-adjoint on 

D(H e ). To see this, note that —ieS/x is bounded with respect to e 2 Ax with relative bound 
and that for ip £ D{H e ) 

||( £ 2 A X ®1)V|| < ci (}\(s 2 D A ®l) 

< 



ci 
C3 



l(^ 



+ 

1 + 1 ® flb 



+ 



< c 3 + W) , (32) 

where we used that H e Q is bounded from below and that H e \ is bounded. Hence H e — i?| iag 
is infinitesimally operator bounded with respect to H £ , consequently H| iag is self-adjoint on 
D(H £ ) and thus @ holds on D(H e ). 
and @) in (H) give 



P. 



-iH% tie 



-l£ e dla f 



t/e 



-iHH/e 
t/e 



(33) 



e"W pj- (V X P*) P* ■ (eVx ® 1) e 



+ 0(e)|t| 



where we used that the term of order 0(e 2 ) in (31) yields a term of order 0(e)|t| after integration, 
since all other expressions in the integrand are bounded uniformly in time and the domain of 



integration grows like t/e. In (^) and in the following we omit the adjoint term from (29) 
and thus consider the difference of the groups projected on RanP^ only. The argument for the 
difference projected on RanP* goes through analogously by taking adjoints at the appropriate 
places. 

Now e(VxP*) • (eVx <S> 1) is only 0(e) in the norm of C(W ,e (&Ti e ,Ti) and thus, according 
to the naive argument, only 0(l)|i| after integration. As in |l3| and p3| we proceed by writing 
(VxP*) • (eVx (8> 1) as the commutator of a bounded operator B with H £ modulo terms of 
order 0(e). This is in analogy to the proof of the time-adiabatic theorem [|lf| and allows one to 
write the first order part of the integrand in (33) as the time derivative of a bounded operator 
and, as a consequence, to do the integration without losing one order in e. 

In view of ((2^) we define 

(34) 



2m 



B(X) :=—.& d\R x (H e (X)) 2 P, L (X) (V x H e )(X) R X (H C (X)) P*(X) 



j(X) 



An easy calculation shows that 



He, P> 



Pj-(V X P*)P* 



(35) 



By assumption d x H e (X) £ C 2 (R n , C(H e )), j = l,...,n, hence Bj(X) £ C 2 (R n ,£(H e )) and 



thus 



-D 2 A ®1,B 



-e (V X B) ■ (eV x ® 1) + 0(e 2 ) = 0(e) 



(36) 
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in the norm of C(W 1,£ £g> Tt e ,Ti). ( |35| ) and (^) combined yield that 

-Pj-(yxP*)P* + 0(e) 



H £ ,B 

with 0(e) in the norm oi C{W l,£ ®TL e ,TL). Since V xH e € C(T~t), a short calculation shows that 
[H £ , eV x ® 1] = O(e) in ^W 1 " 6 ® H e , W)- Hence we define 

B:=B ■ (eV x 8> 1) 

and obtain 

[H £ , B] = - (VxP*) P, • (eV x ® 1) + 0(e) 
with O(e) in the norm of C(W X ' £ ® H e ,H). Let 

B(s) = e iHEs Be~ iHEs 

then 

B(s) = j Ws [H £ ,B] e~ iHEs . 

as 

Continuing (^), we have 

= iee^WI* dse tH ^ s [H £ , B] e~ lHEs + 0(e)(\t\ + \t\ 2 ) 
J o 

= ee -«W/* J^ S dse iH ^ 8 e- mes (J- B(s)) + 0(e)(\t\ + \t\ 2 ) , (37) 

where 0(e) holds now in the norm of ^(W 1 ' 8 ®^, H). The additional factor of \t\ in pi) comes 
from the fact that 

ll^lU^^e) < c(l + eH) (38) 
for some constant c < oo, i.e. the scaled momentum of the nuclei may grow in time. Using 
||-4ext||oo = C < oo and 

||[( £j D A (8)l),B e ]|| £(w) < Ce, 

(pq) follows from 

||HeVx® l)e- iifSs V|| < ||(eZ?A»l)e-^ £s V|| + ||(eA«t ® 1) e^VH 

< IKeDa 1) VII + || [(sD A ® 1), e"^] V|| + C HVII 

< ||(-»eVx® 1)^|| + Ce|s| W + 2C||^|| 

for ipeW 1 ® n c . 

Finally, continuing (p7|), integration by parts yields 

ee -^W/ £ jf*' d8e iH ^'e- 4He ' (4; B ( S )) +0{e)(\t\ + \t\ 2 ) 
e (Be- iHH / £ -e- iH ^ t/e B) 

pt/e 

+ iee -i H ti^/ £ / tfo e «W [ H £ - H £ di )Be~ lHes + 0(e)(|t| + \t\ 2 ) 
Jo 

0(e)(l + \t\f, (39) 
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where 0(e) holds in the norm of C(W 2 ' £ <g> H e ,H). For the last equality we used that B is 
bounded in C(W 2 ' e ® H e , H) as well as in C(W 2 ' £ 7i e , W 1,£ (g> 7i e ) uniformly with respect to e, 
H e - H £ diag is 0(e) in C(W^ £ ® H e ,H), as we saw in (||) and (|T|), and 

\\^ tH£s \\c(w^ Hc )^ ^a+^Ni) 2 (40) 

for some constant c < do. ( [iCQ follows from arguments similar to those used in the proof of (|38l). 

We are left to prove (|l3|). This follows from exactly the same proof using that £(H £ ) 
commutes with e~ tHEs and that, according to (|32"D, 

||( £ 2 A X ® 1)^)^11 < c 3 (||iP £(^ e )V|| + W) < C4(|f| + 1) W- 
4.2 Locally isolated bands 

To prove Theorem ^| we proceed along the same lines as in the previous section, with the one 
modification that we use Proposition ||| to control the flux out of dA. However, one cannot 
use P* = dX P*(X) to define anymore, because the functions in its range would not 

be in the range of H £ and some smoothing in the cutoff is needed. For i G {0,1,2,3} let 
]lj = fl^ A 4 _i S be approximate characteristic functions according to Definition ||. Then the 

smoothed projections are defined with Pi(X) = li(X) P*(X) as P, = J® dXPi(X). In the 
following it will be used that for i < j we have P%Pj = PjP% = P% and hence (1 — P/)Pj = 
P(1-P,)=0. 

Proposition [| yields 

L-iHH/e _ W e -iHg Q t/s j/j pa = f^-iHH/e _ Pj ^* e -«^ *A P° + 0( e ) . (41) 

We make also use of the fact that the phase space support of the initial wave function lies 
in T and has thus bounded energy with respect to H c \. Let E := sup 2gr H c \(z) < oo, let 

^((-oo,_e+o),q) be a smooth characteristic function on R and let £ := ( l((-oo,£;+a),a)(-^ci( - )) ) 
Then standard results from semiclassical analysis imply the following relations. 

Proposition 7. 

(a) l^ a) =£l^ a) + 0(e); 

(b) e~ lH Bo*/e £ = £ e"* H Bo*/e _|_ 0(e) uniformly for t G I; 

(c) [H £ BQ ,£]=0(e 2 ); 

(d) £ G £(L 2 (R n ),W 2 > £ ). 

In (a)-(c) 0(e) resp. 0(e 2 ) hold in the norm of C(L 2 (W n )). 

Proposition (0) (a), (c) and (d) are direct consequences of the product rule for pseudo- 
differential operators (see, e.g., plL|6||) and (b) is again Egorov's Theorem. 
Using Proposition [7| (a) and (b) we continue (|4l| ) and obtain 

L-iHH/e _ pi W e -i£T| t/ e U ^ pa = L-iHH/e _ p lU * £ ^iH^t/e j/j pa + ( £ ) . ( 42 ) 
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We proceed as in the globally isolated band case and write 
[e' iHH ' e - PiU* £ e- iH x° t/£ u) if 

rt/e 

= -ie- lHH/e / dse iues (H £ P x U*e - P X U* £ H^ )e~ iH ^U Jf 
Jo 

rt/e 

= -ie~ lHH/£ / ds e lHEs (H £ - H £ di ) P x U* £e~ lH *o s u If (43) 
Jo 

- ie- iHH/£ / ds e iHSs (H £ di P X U*£- P x W £ #§ ) e - iH Bo s U P% , (44) 
Jo 



where 



H di ag '■- P$H £ P-z . 



One can now show that ( [43[) is bounded in norm by a constant times e(l + \t\) using exactly the 
same sequence of arguments as in the proof in the previous section. One must only keep track 
of the "hierarchy" of smoothed projections, e.g., instead of (29) one has 



(H £ - H diag )P x = (1 - P 3 



e 2 



A x <8> 1, i> 2 



Pi + 0(e 2 ). 



The adjoint part drops out completely, because this time only the difference on the band, i.e. on 
RanPi, is of interest. Note also that the smoothed projections Pj are bounded operators on the 
respective scaled Sobolev spaces and thus, according to Proposition ^ (d), all estimates hold in 
the norm of C{TL). 

It remains to show that also ( f44| ) is 0{e). First note that, according to Proposition ^ (c), 
commuting £ and Pfjo yields an error of order 0(e 2 ) in the integrand and thus an error of order 
0(e) after integration. For <p S W 2 we compute 

{H^P x U*<j>){X) = l 1 {X)E(X)<f>{X) x (X) + l 1 (X) (y(-iVx+A ext ) 2 ^ (X) X (X) 

+ eli(*)(-tW$ (X)-(- i ( X (X),V xX (X)) Hc ) X (X) 
-ie(Vl x )(X)-(-ieV0)(X) x (X) + 0{e 2 ) . (45) 

On the other hand, again for cj) £ W 2 , 

{P x WH^ Q <f>){X) = t l {X)E{X)4>{X)x{X)+t l {X)(^{-iV x + A c ^) 2 c^(X)x(X) 

+ e MX) (-ieV$ (X) • A gco (X) X (X) + 0(e 2 ) . (46) 

Hence 

H £ diag Pi WS - Pi W H{% £ = -eW (Vli) • eV x £ + 0{e 2 ) 
Thus the norm of (H) is, up to an error of order 0(e), bounded by the norm of 

rt/e 

eU* / ds{Vl 1 )-eV x £e-'' H vo s UP 1 ? . (47) 
J o 



(Vli) • eV x £ is a bounded operator and we can apply Proposition |3] in the integrand of (|4 
once more, this time however with the smoothed projection Pq, and obtain 

rt/e 

(47) = eW / ds(Vli) -eVx£%e- lH vo s i{p« + 0(e) = O(e). (48) 
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The last equality in (48) follows from the fact that [eVx £, Ho] = 0(s) and that (VHi) and llo 
are disjointly supported. 

Proof of Proposition |^. For the following calculations we continue %(•) £ C^(A,Ti. e ) ar- 
bitrarily to a function x(') £ C-g°(R n ,H e ) by possibly modifying it on A \ (A — 5/2). For eft 
in a dense subset of L 2 (A — 5) and X £ A — 5/2, by making the substitutions k = ek and 
Y = (Y — X)/e and using Taylor expansion with rest, we have: 

((a w > £ ® 1) X ) (*) = (27r)-» | dY d£; a (^y^, £*) e^*"^* <t>(Y) X {Y) 
= (27r)- n J dYa® (x + -Y, -f) 4>(X + eY) X (X) 

+ e(27r)- n J dYaW (x + |y, -f) <f>{X + sY)Y ■ (V xX )(f(X,eY)) 



= (U*a w > e U(f>x) (X) + R £ . 
From ( |49|) we conclude that 

||(1a-5/2(0®1) (« W ' £ ®1 - W*a W ' £ W)PA- 5 || < P £ ||. 

Since 

||(l-ll A _i /2 (-)®l) (a w ' £ ®l -^a W ' £ W)P A _ 5 || 
= ||(l-l A -«/a(-)®l) (a w ' £ ®l - W*a w ' e W)(ll A _ 5 (-)0l)P A - 5 || =0(e n ) 

for arbitrary n, Proposition ||] follows by showing that i? e is of order e: 



(49) 



(50) 



|i? £ 



< e (2tt)" 



dY 



?( 2 ) 



+ eY)Y-(V xX )(f(;eY)) 



< e(27r)- n sup ||(Vxx)WII«. 

xeR" 



dY 



?( 2 ) 



e y) 



L 2 (R") 



< eC 
= eC 



L2{M „) / dY sup \Y\\a^(X,Y)\ 
J xeR n 



H 
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